Superconducting transition in disordered granular superconductor in magnetic fields 
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Motivated by a recent argument that the superconducting (SC) transition field of three- 
dimensional (3D) disordered superconductors with granular structure in a nonzero magnetic field 
should lie above H C 2(0) in low T limit, the glass transition (or, in 2D, crossover) curve H g (T) of 
disordered quantum Josephson junction arrays is examined by incorporating SC fluctuations. It is 
found that the glass transition or crossover in the granular materials can be described on the same 
footing as the vortex-glass (VG) transition in amorphous-like (i.e., nongranular) materials. In most 
of 3D granular systems, the vanishing of resistivity upon cooling should occur even above H C 2(0), 
while the corresponding sharp drop of the resistivity in 2D case may appear only below H C 2 as a 
result of an enhanced quantum fluctuation. 



PACS numbers: 74.40.+k, 74.81. Dd, 75.10.Nr 



I. INTRODUCTION 



Throughout extensive studies on phase diagrams of 



nongranular system; 



.1.2.3 



it is understood at present 



that the superconducting transition, characterized by 
the vanishing of resistivity, in homogeneously disordered 
(amorphous-like) type II superconductors under nonzero 
fields occurs as a glass transition. As far as the static 
disorder is point-like, the resulting glass transition curve 
H g (T) at nonzero temperatures (T > 0) does not de- 
viate much from the melting transition line of a clean 
vortex lattice. For amorphous-like (nongranular) three- 
dimensional (3D) systems, the H g (T) curve, determined 
resistively, approaches a field near iJ c2 (0)^ in low T 
limit, even if it apparently approaches a field range below 
7J c2 (0)ii upon cooling as the reduced temperatures T/T c o 
is not low enough. On the other hand, It is believed that 
the glass transition in 2D disordered case occurs only 1,6 
at Hsi — H g (T — 0) below -ff c2 (0) and corresponds to 
a field-tuned superconductor-insulator (SI) transition at 
T = oSiS 

However, it is unclear whether this picture also holds in 
granular systems or not. By neglecting effects of vortex 
pinning and superconducting (SC) fluctuations, a glass 
phase peculiar to disordered granular superconductors 
was obtained^ as a phase lying above H C 2(T) definecUS, 
at longer scales than the intergrain spacing (see Fig.l). 
Hereafter, this glass phase, which may appear even in 
H = separately from the Meissner state, will be called 
as the phase glass (PG). However, the fate of PG is not 
clear once SC fluctuation and vortex pinning effects are 
taken into account to describe real systems. It was argued 
recently within a standard model with dissipative quan- 
tum phase dynamics that H g (T) of 3D disordered gran- 
ular systems near T = should lie far above iJ c2 (0)ii*i£ 
and that, even in 2D granular systems with no genuine 
glass transition- in T > 0, the corresponding crossover 
line H g (T) defined from a sharp drop of the resistance 
should show, upon cooling, a divergent (upward) low T 
behavior as if it terminates at a field H g (Q) higher than 
if c2 (0)ii in low T limit (see Fig. 2). This argument may 
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FIG. 1: Mean-fielA H-T phase diagram of a disordered gran- 
ular superconductor. 



normal 




FIG. 2: H-T phase diagram, proposed phenomenologically— , 
of a disordered granular superconductor. In 3D case, the re- 
sistance vanishes on H g (T) (solid curve). 



be consistent with the presence of PG in the mean field 
phase diagram Fig.l if the PG is superconducting. How- 
ever, if so, it is unclear how the SC (glass) phase and 
the portion of H g (T) in H < iJ c2 in Fig. 2 are described. 
Further, the argument in Refiii for 2D case is incompat- 
ible with the field-tuned S-I transition behavior which is 
believed to occur below H c2 (0) even in granular SC thin 
films—. 

In this paper, the glass transition curve H g (T) of disor- 
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dered granular superconductors is examined in the mean 
field approximation but by including SC fluctuations. We 
show that H g (T) in the granular case is obtained in for- 
mally the same manner as the vortex-glass (VG) transi- 
tion curve in the amorphous-like case and hence, is the 
SC transition curve in 3D systems at which the resisi- 
tivity vanishes^* at least in type II limit. Further, we 
find that, even if the mean field PG phase is absent even 
at T = 0, a phase diagram of the type of Fig. 2 with 
H g (0) > H C 2(0) is generically obtained in 3D, while the 
situation in which H g (Q) < H C 2(0) is easily reached in 
2D, reflecting an enhanced quantum fluctuation. We ar- 
gue that the main origin of H g (0) lying above H c2 (0) is 
an enhancement of vortex pinning due to the SC fluctu- 
ation. 

In sec. II, a model of disordered and dissipative 
Josephson-junction array is introduced and rewritten 
into an effective action, and H g (T) curves in 2D and 
3D cases are derived in sec. III. For comparison, the VG 
transition curve in 3D amorphous case is given in the 
same framework as in IV, and related discussions and a 
summary are given in sec.V. 



II. MODEL 



We start from the hamiltonian 



d 



3 V '"'j <i,j> 

describing a Josephson junction array with a charging 
energy 2a on each grain, where the pair of indices, < i, 
j >, denotes a nerest-neighbor pair of sites (i.e., grains), 
and 6j is a phase operator on the j-th grain. Below, the 
model will be extended to a more general one including 
effects of possible dissipation on each grain and of elec- 
tromagnetic fields. The most straightforward method of 
performing this is to express the model into the cor- 
responding quantum action 



5 = 5, 



J ° <i,j> 



J, 



exp[i(^(r)-^(r) 



(2) 



e*5A i -j{T))) + c.c. 



dr / dr' 



[l-co S (^(r)-^(r')); 



^-^ 7t(t — T 1 ) 

j ' 



in the unit % = c = 1, where (3 = 1/T, e* is the Cooper- 
pair charge, A cxt .i-j denotes the line-integral of an ex- 
ternal gauge field over the bond i — j, SAi^j(r) is the 
corresponding gauge disturbance introduced for obtain- 
ing the conductivity in the i — j direction, and 



1 ( d0i{r) 
dr 



(3) 



is the action corresponding to the charging energy, i.e., 
the first term of eq.(l). 

Note that the dissipative (last) term of eq.(J2J) is ex- 
pressed as 



3 <" 

where is the Fourier transform of 

<Z> j (r)=ex 1? (i6 j (T)). 



(4) 



(5) 



That is, eq.JU is nothing but the familiar dissipative 
term, written in the phase-only approximation, in the 
time-dependent Ginzburg-Landau model. 

The quenched disorder in the system is incorporated 
into a randomness of Jij = J*^ with a nonzero real mean 

Jo, i.e., Jij — Jo > 0, and a Gaussian distribution 
{Jij — Jo)(Jji — Jo) = J 2 - These relations may be re- 
garded as being due to a random gauge field ay defined 
by Jij — Jo oc exp(iay ). The free energy F = — /3 _1 lnZ 
will be expressed in terms of the replica trick as F = 
—(3~ 1 (Z n — l)/n in n — > +0 limit. The averaged repli- 
cated partition function Z n is given byi£ 



Z n = Z% < exp(-5/ - S g ) > , 



(6) 



where Zq is the partition function of 5o, < >o denotes 
the ensemble average on Ei<a<n^o(^°')' an d 

n .0 

S f = "E E / dTj cos(e*6A i _ j (T) + e*A cx , l _ J (7) 

a=l <i,j> ^° 

- ef\r) + ef\r)) + J dn J 
x [l-cos^^n)-^^))], 

S g = -\f dr f dr'Y^Y, -PanWAi-jfj) 

a,b <i,3> 

- SAi-iV)) + <tV) - 9f\r') - 9f\r) + ef\r')). 



Before proceeding further, 5/ will be rewritten in the 
fori'" 11. 



s { = const.-dJor'EE^H)^ 1 -^^! ( 8 ) 

1 



2d 



D 4 .Dn$| a) H 



const. - d/o/T 1 E I 1 + " M/(2dJo)] _1 ^($f>)) 



(l + -D l D; (a,) 



(9) 



for the cubic or square lattice in e?-dimension, where 
$,( T ) = (3' 1 XL $i{u)e~ iuir . Equation (9) is valid up to 
the lowest order in v\u>\/ J$ and the laplacian Dj-D?/(2<i), 
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and Di is the gauge-invariant gradient on the latticed 
accompanied by the gauge field A cx + <5A(r). Then, 
by introducing the conventional SC order parameter 
ipi (r) and the glass order parameter q^ (ri, t 2 ) = 
(g^(r 2 ,ri))*, Z" becomes&iSiia 

7n r 

= = v^v{^yv q ^eM-s c ^,q)l (10) 

where 

S cS ty,q) = fdrJ dr 2 J2J2^r^ ab) ^ T ^ C 11 ) 



we finally obtain the following effective Landau action 



<9t 2 



(14) 



x Q( aa )(x;r 1 ,r 2 ) 



a, 6 i 



u) i a 



2d J 



(«)|2 



_ »y dndradra £ Q< a6 >(x; n , r 2 ) 
Q (&c) (x;r 2 ,T3)Q^(x;r 3 ,r 1 ) + ^ / dr ^(Q( aa )(x; r, r)) 2 
J2j dT iJ dr 2 \(-iV - e*(SA( n ) - SA(t 2 ))) 

tS e ff, 



in 



- ^ln^r T exp(^yrfr^$W(r)(l 



4da 2 

x Q( afc )(x;r 1 ,r 2 )| 2 
where 



i 5 c ff — a 



/ dsJW/r^^Miv^i 2 ) (15) 

** L „ \ ..... 



2d 



X ($f ) (T 2 ))*ll + 



dr 



^,o|^ (a) (r)| 2 +c^ 



9-0 



(a) 



9r 



d.d; 

2d 



1/2 



q\ (t 2 ,ti) + c.c. 



to'K-iV - e*A ox - e*<yA(T))^w(T)|" 

2 



where 



2a \ a 



4Jo 

a 



(12) 



l^ (Q) (x,r)| 4 
" dn y dr 2 (^ Q )(x,r 1 ))*Q( Qb )(x;r 1 ,r 2 ) 



and Di denotes the gauge-invariant gradient on the lat- 
tice accompanied by the gauge field <5A(ti) — SA(t 2 ). 
Performing the cumulant expansion in powers of 
and ip^ a ' in the logarithmic term, various terms such as 



dnl dr 2 <T r $( a )(r 1 )(<i>( Q )(r 2 ))* > g (aa) (T 1 ,r 2 )(13) 

IJ dn J dr 2 J dr 3 y dr 4 <T t $W(t 1 )($( q )(t 3 ))* 
X $W( T4 )($«>) ( T2 ))* >og {6o) (7a,ri)g^(7a,r4), 



arise in the resulting Landau action S c g appropriate to 
the ensuing analysis. The average < >o is carried out 
by using So or its soft-spin version^. For instance, < 

T T ^ a) (n)(^ a) (r 2 ))* > becomes exp(-a|n - r 2 |/2) in 
the low T limit. Below, it will be replaced by its local 
limit 4(5(r)/a anywhere except in the lowest order term 
in q( ab K Further, the T-dependence will be taken into 
account just in the I-0I 2 term because, at least, one of such 
T-dependences is necessary in order to keep a reasonable 
mean field H c2 (T) line for the ^-field. It will be clear 
that these simplifications are not essential to the present 
purpose of addressing the low T phase diagram. 

Next, let us write q^ (t\ , t 2 ) , by following Read et 
al4£, as Q < - ab) (ri,T 2 ) — C8 ai bS(ri - r 2 ) in order to delete 
the term J dr\ J dT 2 \Q( ab \ri , r 2 )| 2 . By representing spa- 
tial coordinates in terms of the continuous coordinates x, 



X ^ (b) (x,T 2 ) 



Here, we have introduced a short length cut-off a which 
corresponds to the intergrain spacing. We assume that 
a is much longer than the coherence length of the host 
material forming the grains, and hence that the averaged 
H c2 (§) of the granular system is lower than the micro- 
scopic H c2 (0) of the host material forming each grain 10 . 
The H c2 (0) mentioned in sec. I is nothing but this aver- 
aged H c2 (0). This is consistent with the assumption in 
choosing the phase-only model (1) as a starting model 
that the amplitude of the pair-field in each grain be ro- 
bust. 

We note that, although the dissipative term in eq.Q 
is reflected only in the term quadratic in ip of the effec- 
tive action, the dynamics of the glass fluctuation SQ^ ab ^ 
also becomes dissipative through the coupling (w^-) term 
between t/j and Q^ ab ^ after integrating over the SC (ip-) 
fluctuations. 

Using the soft-spin version of the zero-dimensional ac- 
tion Sq, i.e., a Ginzburg-Landau action corresponding to 
SqI^, we find the coefficients to be given by 



t 



T' 



(16) 
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1 



a 
2J 



>>,o 
t ~ 

T ~~ 



h 2 



— I I — I +1 

2J 2 
a 2 

JqT 
2a 2 



JL _ A 

A 

2a 3 ' 

tv 



■h 



4a V V 2 J 



- 1 



8d 2 J ' 

2 ^0 



a 



4da : 



Jot 



The coefficient u# denotes the renormalized four-point 
vertex in the soft-spin version of Sq. In low dimen- 
sional cases with d < 2, the renormalization of the fluc- 
tuation in a classical Ginzburg-Landau action is well 
approximated 2 ^ by the Hartree approximation in which 
Ur = 0. Based on this fact, the vertex un/a of the quan- 
tum zero-dimensional action So may be assumed to take 
a (dimcnsionlcss) number much less than unity. 

Below, the Fourier transform of the glass field Q( ab )(x) 
is defined, by following Ref>i&, as 

Q( ab )(n,r 2 ;x) = ? (ai) + ^ 5 « e " iu(T1 " ,s) ^( 17 ) 



— 1CJ1 T1—IUJ2T2 



where the replica symmetric form 



>*>) _ 



= <?(1 - S a , 



is assumed for q ab , because we do not study here the glass 
phase below H g . Further, we focus hereafter on the situa- 
tion that H g is approached from the higher temperatures 
at which q = 0, and q = /3 -1 Z?o^- The replica-diagonal 
component D u is the Fourier transform of the (imagi- 
nary) time correlation between two "spins" (Re<F,-, Im&j) 
and hence, is nonvanishing even above H g (T). It is de- 
termined by the variational equation = n~ 1 dZ n / dD u]1 
or 

kdI = K- 1 {LU 2 +r)+uf3- 1 Y,D L ,-w lf ,(3- 1 (^*J^ UJ (x.)) s , 

(19) 

where ( ) s denotes the ensemble and space averages. The 
physically meaningful solution of eq. 1)19(1 is 



where 



UK 



A, = -K- 1 Vr^+tu 2 , (20) 
^ K /3- 1 (C(x)Vo,(x)) s .(21) 



-1 \ 



The minus sign of eq. 1(20(1 is chosen so that a physically 
correct " spin" correlation along the (imaginary) time di- 
rection is recovered^. We note that the last term of 
eq. (|21|l is nonvanishing in T — » limit. 



III. GLASS TRANSITION LINE 

In this section, we will determine the glass transition 
field H g (T) of the granular system described by eas. 1(14(1 
and l(15|). It will be seen below that, upon cooling, the 
glass transition is described as a vortex-glass orderingii^ 
induced by the coupling between the SC fluctuation ip 
and the glass fluctuation 5Q^ ab \ 

We will focus on the (highest) instability temperature 
at which the glass fluctuation 5Q^ ab ^ becomes critical at 
the Gaussian level when the SC fluctuation is fully incor- 
porated. It corresponds to a mean field glass transition in 
analogy to the normal to Meissner mean field transition 
following from the BCS theory (Note that the quasipar- 
ticle in the normal state in the latter corresponds to the 
SC fluctuation in the former). For simplicity, we shall 
identify this mean field glass transition line with H g (T). 
Further, since our main purpose here is to give a correct 
answer on H g (T) in high H and low T portion of the H-T 
phase diagram, we will use the lowest Landau level (LLL) 
approximation for the ip modes. For the LLL modes, the 
operation (— iV^ — e*A cxt ) 2 is replaced by |e*|if , and the 
number of field-induced vortices may be expressed as the 
total magnetic flux HS multiplied by |e*|/27r, where V~S 
is the linear system size. 

Under the nonzero D u given by ea. 1(20(1 . the glass 
fluctuation SQ^ ab ^ obeys the following effective action 
5S e e = SSq + SS^ up to the quadratic order in 8Q^ ab \ 
where 



(18) SS Q 



0- 



EE 



a, b u)\.uJ2 

x -V(5Q (ba) (-w 2 ,-wi) 

.(ah) 



J"V V<5Q (a6) (wi,w 2 ) (22) 



t 



(D Ul + D^SQW (wi, wsJSQ^ (-W2, -wi) 



a 



0-1,-1 -d 
p t a 



E 



J ta 2 



J Q a- 2 tD u 



+ c^ 2 + ^MM a) | 2 + ^HH V - e*A cx )^| 2 



E 



tu R J§ 
2(3 2 



2 -u. 

a 



u)\ -\-u)2 ,cJ3+a;4 



The uSQ^SQ^ term (see eq-lJUJ) was dropped from 
5Sq. In fact, this term actsi^ as an interaction term be- 
tween SQs and hence, may be neglected at the present 
stage of focusing on the noninteracting (Gaussian) 5Q- 
fluctuation. In addition, since ua< yfr for any J/a val- 
ues of our interest, this term is quantitatively negligible 
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even in obtaining (see ea . (|2*TJ) 1 . where denotes 
D u in G?-dimension. Hence, let us drop this small term 
consistently hereafter. Then, the T and H dependences 

in d[j ^ arises primarily from the SC fluctuation. 

First, let us give the renormalized -0-fluctuation in or- 
der to reasonably describe the field range lower than 
H C 2(T). Hereafter, the |-0| 4 term will be treated, as in 
the nongranular casa2£, in the Hartree approximation. 
The Hartree approximation in LLL may be invoked as 
an infinite-range limit of a nonlocal Ginzburg-Landau 
modelSi. Although the use of the nonlocal model makes 
the details of vortex positional ordering obscure 2 ^ the 
description of the glass ordering in strongly-disordered 
superconductors is not essentially affected by this proce- 
dure. Then, the ^-propagator in LLL is defined in 2D 
[d = 2) case by 



(2) 



dial 



where 



h& = 



\e*\Ha 2 
2tt 



(2)' 



Jo V 1 t * 



2^( 2 ) 



(23) 



(24) 



(25) 



(26) 



a- l D {2 ) + £< 2 > 



and 



K.t 



(2) 



1/2 



(27) 



In 3D, the corresponding expressions are given by replac- 
ing G ( 2 with afdk z /(2Tr)G < § a (k z ), where 



Jo(nl 3) +fc 2 a 2 /12)' 



(28) 



Jo V t t t 



2V3^ 3 > 
36 



(29) 



D {3) = -n- 1 \r + uj 



nt 



1/2 



and 



h& = 3 1 / 



1/2 \e*\Ha 2 



2n 



(30) 



(31) 



The self energies T,^ and T,^ due to the interaction 
(ur-) term are given by 



a * — ' 



and 



E (3) =32 ^ (/3a) -^(3)^ (/ ,(3) r l 

a L — ' 



/2 



(32) 



(33) 



It is easily verified that, in J -C a, eas. (|26|l and H29|) 
reduce to their results in the J = case. 

Although solving exactly these set of equations for 
each d is in general difficult, it can be performed just 
at the (mean- field) glass transition line H g (T) . To define 
H g (T), let us first rewrite ea. l|22[l into an effective action 
SS c s,q consisting only of SQ. Within the Hartree ap- 
proximation for the ^-fluctuation, <5<S e ff,Q in 3D becomes 



5S, 



Off,! 



(34) 



a, b oj\,u)2 



where v^ ± = exp(— (kf + fc 2 )/(2/?/ 2 ))). Then, by focusing 
on the term with uij = and k = in eq.lJSi}, H g (T) is 
defined as 



2t" 1 J D 



(3) 
;=0 



ft( 3 > 

2na 2 



(3)\-3/2 



00 



in 3D, while 



- 2t- 1 D {2} 



na 



(35) 



(36) 



in 2D, where n Q d) = mL=o 



To rewrite these equations 
explicitly, /Lt^ at £f s (T) will be expressed as /i^ + 
5fj,jy. If using ea. (|26|l or l|29|l . it is not difficult to obtain 

as a funct: 
u>\, it becomes 



as a function of /Uq , and, up to the lowest order in 



(3) 



8(i 



in 3D, and 



(3) J 

fi d^a 



3 + 64(4 3) ) 5/2 /^ (3) tJ ° 



1/2 



(2) 

Mo 



8(m 2) ) 3 / /i(2) tJ ° 



1/2 



in 2D, respectively. The expressions on obtained 
above imply that, to the lowest order in the Matsubara 
frequency, the glass fluctuation propagator just at the 
(mean field) glass transition takes the form 

(<JQ£ a VV)^W,-")> * (fc 2 + |c| 1 / 2 + | w '|V2 ) -i 

(39) 



(37) 



(38) 
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FIG. 3: Examples of #ig (T)-lines in 3D case following from 
the present theory. The H C 2(T) line is given by h — 1 — 
T/T c o. See the text regarding the parameter values used for 
calculations. 




FIG. 4: Examples of Hi-g J (T) in 2D (sharp crossover lines) 
comparable with those in Fig.l. See the text regarding the 
parameter values used for calculations. 



after being rescaled spatially. 

Now, by applying the above expressions of Sfj,^ to the 
selfconsistent equations on jj^ 



(2) at 
Mo = T7?>,o 
Jot 



and 



(3) a 

Mo = -77^,0 + 
Jot 



2nh<- 2 



36 



a D,., =l 



+ E (S ', (40) 



-17T (3) 
a D,.,- n 



<3) 



(41) 



we obtain the coupled equations, eas. (|26[l . i|27|) . (|32[1 . 
(ESI, (EH, and J30I) in 2D and eqs.®, ED), CSJl, CH, 
(|37|l . and (|41() in 3D. The resulting iJ-T relation for each 
d is nothing but the H g (T) line. 

Typical examples of computed H g (T) lines are shown 
in Fig.3 for 3D and Fig.4 for 2D, where h = H/H c2 (0), 
and T c o is the zero field transition temperature. We 
have commonly used the parameter values, 2ird^Jo/t = 
27TI//32 = 10 and un/a = 1 x 10~ 3 , and have changed 
jo = Jo/ot and j = J /a. In Fig.3, the pairs of param- 
eters (j, j ) are (0.05, 0.7) (left curve) and (0.28, 0.7) 
(right), while in Fig.4 they were chosen as (0.003, 0.3) 
(left curve), (0.05, 0.7) (center), and (0.05, 0.3) (right). 
The dashed curve in 3D is the thermal melting line in 
clean limit. 

In 3D case, the resistivity should vanish at H g (T) , be- 
cause, as is explained in sec. IV, the obtained H g (T) is 
essentially the same as the VG transition fielcUii^ in the 
amorphous case. In contrast, H g {T > 0) in 2D will be 
regarded as a crossover line along which the resistivity 
remarkably dropsii, while H g (0) may be identified with 
the SI transition field Hsi at T = (see sec. I). 

In 3D, the low T limit H g (T — > 0) of the glass transi- 
tion field lies, for most of parameter values we have ex- 
amined, in h > 1, i.e., above the mean- field H C 2(0). On 



the other hand, all values of the disorder strength j used 
in the figures satisfy |Z? w= o| > 0. This implies assuming 
that, even at T = 0, the PG order in H > H c2 (0) does 
not occur. That is, Fig.3 implies that H g can lie above 
H C 2(0) more easily than the prediction from the mean 
field analysis** The presence of H g (T) above H c2 (T) at 
low T is due primarily to the coupling, appearing through 

,D^2o> between the glass field and the SC fluctuation: For 
brevity, let us imagine that the SC fluctuation term in 

(D^1 Q ) 2 (i.e., its second term) will be relatively small in 
magnitude. Then, from ea. (|35|l or (|36(l . /i^ at H g (T) 
approximately scales like (ft( d )) 2 /( 6-ci ). Thus, the sec- 
ond term, - /i( d )/( Mo rf) )( 4 - d )/ 2 , in (D^l ) 2 grows like 
(/j( d )) 2 /(6-d) w j^h increasing field h^ d \ As is explained 
in sec.V, this implies that the strength of vortex pinning 
is enhanced by the SC fluctuation with increasing field, 
although the relation ~ (/ l ( d )) 2 /(6-d) itself implies 
a reduction of the SC fluctuation with increasing field 
above H c2 . Such an enhancement of pinning due to the 
SC fluctuation is peculiar to the granular superconduc- 
tors at low T and in high fields. In fact, at high T and 

in weak or intermediate fields, and hence, |-D^2 I 
rather increases with increasing temperature. 

Further, it can be seen from the coupled equations 
leading to the figures that, in the present quantum model 
in which the main bare energy scale is not the Joseph- 
son coupling Jo but the charging energy a, the main 
Jo-dependence appears in d^a/(Jot) cx j^ 2 , and hence 
that a larger j leads to an enhancement of quantum 
fluctuation 23 ' 24 . For this reason, the (0.05, 0.7) curve 
lies below the (0.05, 0.3) curve in Fig.4. 
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IV. REVIEW OF RESULTS IN 
NONGRANULAR CASE 

For the purpose of understanding the content of results 
in sec. Ill better, it is useful to compare the results in 
granular case with those in nongranular cas o 14 i 25 . Here, 
we will sketch the corresponding analysis for obtaining 
the VG transition curve H g (T) of the amorphous-like 
materials. Within LLL, the familiar GL action derived 
microscopically takes the form 



vanishes in the limit of vanishing k and ujj . In the present 
case, the VG transition line is given by 



d 6 r / dr 



(*(r,T))*u(Q;r)tf(r,iQl2) 



+ £o 2 |3 2 *(r,r)| 2 + -|vI/(r,T)| 4 



j_ f |^(r,T)-vE(r,r')| 2 
2tt J (t-t') 2 



Sgl.q - J dT J dT 'Y,J d 

a,b 



1 



2Ar 



SQ^^t.t') (44) 



x SQ iba) (r; r', r) - (* (a) {t))*SQ^ (t, r')# (b) (t') 

after integrating over SQ. If restarting from the action 
>5gl[^^]+<5gl,0) we can follow a similar route for 
obtaining the glass transition line to that in sec. III. First, 
the propagator of renormalized ^-fluctuation in LLL is 
specified by the Matsubara frequency u> and, in 3D case, 
the wave number k in the direction of the applied field, 
and given in the form (// + j\tv\ + ^fc 2 ) -1 in the Hartree 
approximation 2 ^. For weak disorder, the parameter /i 
satisfies 



Inh + bh^ 3 T 



1 



vm + tM' 



(45) 



which corresponds to ea. H41l) in sec. III. The VG transi- 
tion field can be defined, just as in ea. i|34|) . as a critical 
point of the glass fluctuation SQ at which the inverse 
of the propagator (<5Q( ab )(k; iv 2 )SQ^(k; -u) 2 , -wi)) 



1 



Ah 

2 M 3 / 2 



(46) 



where u — Inh, h = H/H c2 (0), (u(r) — u)(u(r') —u) = 
Ao6^ 3 '(r— r'), b > 0, 7 > 0, and Q is the gauge-invariant 
gradient in directions perpendicular to the applied field 
|| z. The role of the random potential leading to the 
vortex pinning is played by u(r) — u. The replicated 
action, arising after the random-average, takes the form 

n 

S G L,r = £).Sgl[* (o) ] (43) 
a=l 

- ^ J d 3 rY, J dT J rfr'|*( Q )(r,T)| 2 |*W(r,r')| 2 . 

a, b 

Alternatively, the last term of <Sgl,i- may be regarded as 
arising from the action 



which corresponds to ea. (|35|l in the granular case, where 
A = A /(2< 3 ). 

In previous worksiiS, eq.® was obtained as a 
pole of the VG susceptibility directly constructed from 
the action H43|) . Then, the VG susceptibility takes 
the same form as the glass fluctuation propagator 
(5Q( ab )(k;wi,w 2 )(5Q( fcQ )(k;-W2,-wi)) in the limit of 
vanishing k and u>j. That is, if starting the analysis in 
sec. Ill from an effective action composed only of corre- 
sponding to eq. 143|) , the glass transition and resistive be- 
havior near H g (T) for the granular case can be described 
in the same manner as those performed elsewhere 1 * 2 ^. 
Since the continuous vanishing of resistivity at the VG 
transition was explained there based on ea. H43l) . H g (T) 
for the 3D granular case in sec. Ill has to be also the SC 
transition line at which the resistivity continuously van- 
ishes. 

The VG transition field in 3D case and at T — is 
determined by the following expression which is obtained 
from eas.p5)l and 



ln/i 



- 1 £ 3 bu; c h 



(/iA/2)!/3 + ( 7 ^ c + (^A/2) 2 / 3 )V2 

hA^ 2/3 
2 



(47) 



where oj c is a high frequency cut-off. We note that the 
fluctuation-corrected (i.e, renormalized) value of H C 2, 
H c2 , is nonzero in 3D systems at T — and is given 
by eq.(H3> with A = 0. Due to this fact and the A- 
dependence of eq. l|47|l . the 3D VG transition field at 

T = is always higher than H^(0) and approaches 
Hj2 (0) as the disorder A diminishesi* 2 ^. In contrast, in 
2D case, Hffl (0) vanishes, and hence, the VG transition 
field at T = (i.e., i?si) m ay lie below the mean field 
H c2 (0). 



V. SUMMARY AND DISCUSSIONS 

First, let us start from explaining the 3D H g (T) go- 
ing beyond the mean field H c2 upon cooling in sec. Ill 
on the basis of the results in sec. IV for the nongranular 
case. One reason for the if g -growth in higher fields seen 

in 3D case is that H g (0) > i?^(0) in any 3D case, and 

that H^ 2 (0) is nonvanishing and tends to lie near the 
mean field H c2 (Q). For instance, even in the left curve 
in Fig. 3 where T„(H) is significantly lowered, H^'{0) 
lies below but close to H c2 (0). For completeness, we 
show the lowest T behaviors of H g (T) line for extremely 
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1.0 




0.002 T/Ta 



FIG. 5: 3D H g (p) curves at very low t = T/T c0 for the weak 
disorder cases with j = 1 x 1(T 4 , 2 x 1(T 4 , and 2 x 1(T 3 . In 
common, the value jo = 1.5 was used. 



low j values in Fig. 5. It shows that, as mentioned be- 
low ea. (|47|l in the nongranular case, H g (0) approaches 

H ( C 2 ] {0) = 0.978 H c2 (0) as the disorder j diminishes. 

That is, since H£\0) is the lower limit of the 3D SC 
(i.e., glass) transition field at T = 0, H g (T — ► 0) in d = 3 
lies in /i > 1 in most cases (see also Appendix). In con- 

trast, H£'(0) in 2D case is zero, and consequently, if si 
of a system with strong enough quantum fluctuation, as 
the left curve in Fig. 3 shows, can lie in h < 1, i.e., much 
below the mean field if ^(O). This situation corresponds 
to the case in which a field-tuned S-I transition behavior 
is seen through resistivity curves in granular materials 13 . 

More importantly, in the granular case at low T, there 
is a contribution of the SC fluctuation enhancing the vor- 
tex pinning strength with increasing field: As explained 
in sec. Ill, the l.h.s. of each of eas. (|35|l and (|36|l decreases 
with increasing field due to the SC fluctuation. On the 
other hand, by comparing ea. l|35|l with ea. (|46|l . the in- 
verse of the l.h.s. of ea. l|35|l . increasing with field, corre- 
sponds to a strength of vortex pinning inducing the glass 
transition. This unfamiliar fluctuation effect, peculiar to 
granular systems, is quantitatively weakened in 2D case, 
because the quantum fluctuation is stronger in lower di- 
mensions and selfconsistently increases /Zq . 

In the present paper, we have examined the low T be- 
havior of SC glass transition curve H g (T) of granular 
superconductors by applying a theory of quantum spin- 
glass to the context of superconductivity and have shown 
that, in contrast to the situations in amorphous-like ma- 
terials, a situation with H g > H C 2 at low T usually oc- 
curs in 3D granular systems. This is consistent with the 
phenomenological picture2iAi*i£ favoring the presence of 
a superconducting glass in if > ff C 2(0). Next, we have 
shown that, without the PG order, the fluctuation of the 



glass order parameter plays the role of pinning disorder 
inducing a vortex-glass instability at T > in 3D systems 
under nonzero fields. 

Our main messages based on the results in sec. Ill 
and IV is that the SC transition in nonzero fields 
can be described in a single theory, i.e., as a vortex- 
glass transition^, for both the granular and nongranular 
(amorphous-like) materials. Because the granular and 
amorphous-like systems may be continuously connected 
with each other, e.g., by changing a composition of mate- 
rials, such a unified view of two limiting models of disor- 
dered superconductors should be naturally expected. In 
a work with a similar purpose to the present one, Galit- 
ski and Larkinii have argued that even the SC transition 
in amorphous- like materials should be described within 
a model for granular systems. Our result in sec. Ill that 
H C 2(0) < H g (0) < +oo is consistent with the argument 
in Refsiii andi£. However, our result in 2D case, given 
in Fig.4, that H g {T -> 0) < if c2 (0) is different from 
their opinion and rather consistent with experimental 
facts showing the S-I transition behavior. 

The vanishing of resistivity on approaching H g from 
above should imply that the glass phase in if < H g is su- 
perconducting, because the transition at H g is continuous 
in the present case. It is not surprising that the present 
result disagrees with an argument^ based on almost the 
same model that even the 3D VG is a metal, because 
the dissipative term of eq. (J2J) was not taken into account 
correctly in Ref»2£. Without the dissipative term, a finite 
fluctuation conductivity at finite T above H g cannot be 
discussed^ In fact, it is difficult to reconcile the metallic 
response^! in PG phase with the conventional ac Meiss- 
ner response^. Further, if their model with no Ohmic 
dissipationSI is extended to the case with correlated line- 
like disorder creating the so-called Bose glass phased, a 
(Bose) glass phase with finite resistivity seems to be still 
obtained in contrast to experimental facts. Our analysis 
on the conductivity in glass phases including these results 
will be presented separately^ 
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VI. APPENDIX 

Here, we estimate H g (0) in 3D amorphous- like case 
at low T limit by using microscopic parameters in dirty 
limit. In the ordinary dirty limit where Tr -C 1, the 
coefficients in the action (|42H are availableS&iffliSi. If rea- 
sonably choosing jlo c to be a constant of order unity, we 
have uj c ~ t _1 , where r is the elastic scattering time of 
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quasiparticles, and 

where £g is the T = coherence length in clean limit, 
T c is the zero field transition temperature, and h depen- 
dences in the coefficients were neglected by assuming the 
glass transition field H g (0) in T — > limit to stay close to 
H C 2(0). This is justified as far as both buu c /£^ and A< 1. 
If T c t < 1, 6w c /Co > A 2/3 > and consequently, -Hg(O) lies 
below H C 2(Q). By contrast, for larger tT c values of order 
unity, H g (0) might lie rather above H C 2(0). Nevertheless, 
its difference (H g (0) — -ff C 2(0))/-ff C 2(0) is small according 



to 

H g (Q) - H c2 (0) ~ H c2 (0) (J) 1 ^ m - (49) 

Of course, when T c t 3> 1, the rhs in ea. H47|) is negligible, 
and H g (0) becomes^ of the order of the nonvanishing 

h£\o) just below the mean field H C 2(0). 

In the above analysis, effects of Coulomb repulsion 32 
between quasiparticles were neglected. As in 2D 
cas o''9'iffl , it would play a role of reducing H g (0), al- 
though H g (0) does not become lower than H^'(0) in 
3D case. 

In this manner, the statement in Introduction that 
H g (0) of nongranular superconductors lies close to H c2 (0) 
is justified. 
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